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Abstract 

Massive gravity is a modified theory of general relativity. In this paper, we study, using 
a method in which the scale factor changes as a particle in a "potential", all possible cosmic 
evolutions in a ghost-free massive gravity. We find that there exists, in certain circumstances, 
an oscillating universe or a bouncing one. If the universe starts at the oscillating region, it may 
undergo a number of oscillations before it quantum mechanically tunnels to the bounce point and 
then expand forever. But going back to the singularity from the oscillating region is physically 
not allowed. So, the big bang singularity can be successfully resolved. At the same time, we also 
find that there exists a stable Einstein static state in some cases. However, the universe can not 
stay at this stable state past-eternally since it is allowed to quantum mechanically tunnel to a 
big-bang-to-big-crunch region and end with a big crunch. Thus, a stable Einstein static state 
universe can not be used to avoid the big bang singularity in massive gravity. 
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I. INTRODUCTION 



The current accelerated cosmic expansion was discovered firstly from the Type la 
supernovae data [f, 2] and further confirmed by many other observations, including the 
cosmic microwave background radiation [3], the large scale structure [4, 5], and so on. This 
discovery broke our common belief that the universe should be undergoing a decelerated 
expansion. A possible explanation for it among other such as the cosmological constant 
and the dynamical dark energy is that the theory of general relativity is no longer valid on 
the cosmological scale and needs to be modified. As a result, many modified gravity the- 
ories have been proposed to explain the present accelerated cosmic expansion without the 
need of the mysterious dark energy. Among them, the Dvali-Gabadadze-Porrati (DGP) 
model [6] is a very interesting one since it not only admits a self-accelerating solution 
with only pressureless matter, but also allows the graviton to have a small mass on the 
cosmological scale. 

Actually, about eighty years ago, Fierz and Pauli [7] first tried to build a theory of 
massive gravity. However, the linear Fierz-Pauli theory can not recover the linearized 
Einstein gravity in the limit of zero graviton mass and can not pass the solar system 
tests due to the van Dam-Veltman-Zakharov (vDVZ) discontinuity [8]. With the help of 
Vainshtein mechanism, the introduction of nonlinear interactions can cure this discontinu- 
ity [9]; unfortunately, at the same time, the nonlinear terms also yield the Boulware-Deser 
(BD) ghost since more than two time derivatives are contained in them [10, 11]. In order 
to construct a consistent theory, nonlinear terms should be tuned to remove order by 
order the negative energy state in the spectrum [10]. Recently, a ghost-free nonlinear the- 
ory of massive gravity was constructed successfully by de Rham, Gabadadze and Tolley 
(dRGT) [12]. (See, however, [13] for the causality issue of the theory.) It has been found 
that the dRGT theory can accommodate cosmological solutions with self acceleration [14] 
and the observational constraint on it has been discussed in [15]. In addition, the Einstein 
static state (ES) universe in this massive gravity theory was analyzed in [16] and it was 
found that there exist stable ES solutions to avoid the big bang singularity problem. Let 
us note that the Hawking-Moss instanton in massive gravity has been studied in Ref. [17]. 
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If one writes the Friedmann equation in a form such that the evolution of the cosmic 
scale factor can be treated as that of a particle in a potential, then it is possible to classify 
all cosmic evolution types as has been successfully done in the Horava-Lifshitz gravity [18] 
and the DGP braneworld scenario [19]. In this paper, we plan to study all possible cosmic 
evolutions in the massive gravity with this method. The paper is organized as follows. 
In Sec. II, we give the Friedmann equation of the massive gravity and define all possible 
cosmic evolution types. In Sec III, we derive the conditions for the ES solution. We then 
classify all the cosmic types and give the conditions for them in Sec IV, V and VI, and 
conclude in Sec VII. 



II. THE FRIEDMANN EQUATION 

We consider the theory of massive gravity proposed in [20]. The action has the form 
S = [ ( - \r + m 2 C J ^jd 4 x + S m , (1) 
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where G is the Newton gravitational constant, R is the Ricci scalar and hm/c is the 
graviton mass. In the present paper, we let S m describe the ordinary matter plus a 
possible cosmological constant generated by vacuum energy. C is the nonlinear higher 
derivative terms for the massive graviton and it is defined as 

r l/o2 qA q B\ i C 3 ABCQqMqNqP i C 4 ABCQ nM qN qP qQ / \ 

L=-[i3 - b B b A ) + — e M NPQt b A b B b c + —e M NPQ£ b A b B b c b D , {I) 
where S = S A , c 3 and c 4 are two constants, €mnpq is the Levi-Civita tensor density and 

si = 6i-yi. (3) 

Here j B is defined by 

lh C A = 9 A °fAB (4) 

with f AB being a symmetric tensor field. 

The Robertson- Walker (RW) metric for a spatially homogeneous and isotropic universe 
can be written as 

ds 2 = dt 2 - a 2 (t) ( + r 2 d 2 n) , (5) 
V 1 — kr 2 / 



where a is the cosmic scale factor, t is the cosmic time and k = 0, ±1 is the constant 
curvature of three dimensional space. In massive gravity, Chamseddine and Volkov show 
that there exist cosmological solutions where the effect of the graviton mass is equivalent 
to introducing to the Friedmann equation a matter source that can consist of several 
different matter types besides a cosmological constant term [21] 

tt2 . k m2 ( a j_ a i Q ^3-3c 3 -c 4 2 c 4 + 2c 3 -l 
H + — = — 4c 3 + c 4 - 6 + 3C h 36 
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where H = | is the Hubble parameter, C is an integration constant and p is the energy 
density of ordinary matter plus vacuum energy. Here, besides the cosmological term, three 
additional terms in the Friedmann equation which decay as 1/a, I /a 2 and 1/a 3 can be 
viewed as quintessence, gas of cosmic strings, and non-relativistic cold matter respectively. 
Thus, the massive gravity can explain the present accelerated cosmic expansion. At the 
same time, it may also play an important role in the very early universe (when a is very 
small). In the present paper, we plan to study the all the possible cosmic evolutions 
in massive gravity and whether the big bang singularity can be avoided. Note that, 
by employing the familiar canonical quantization procedure in massive gravity for an 
open cosmic background, Vakili and Khosravi found that the big bang singularity can be 
avoided through a bounce [22]. For simplicity, we only consider a spatially flat universe 
k = and a positive constant C. In addition, since we are interested in the very early 
universe, the vacuum energy is assumed to be the only cosmic energy component and then 
p is a positive constant. Defining A = and rescaling a/C — > a, we can re-express the 
Friedmann equation (Eq. (6)) as 



m 2 ( 3 - 3c 3 - c 4 c 4 + 2c 3 - 1 c 3 + c 4 

— 4c 3 + c 4 - 6 + 3A + 3 V 3 — 

3 \ a a z a 6 



Let us now write the above Friedmann equation into the following form 

a 2 + V(a) = 0, (8) 
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where 



V(a) = -m 2 -(4c 3 + c 4 - 6 + 3A)a 2 + (3 - 3c 3 - c 4 )a + 

O 

N 1/ , 1 _ 

(C4 + 2C3-1)--(C3 + C 4 )- 



(9) 



Thus the evolution of the scale factor a can be considered as that of a particle moving 
in a "potential" V . Obviously, this "potential" must satisfy the condition V(a) < 0, and 
this gives the possible ranges of a as the universe evolves. Since the values of A, C3 and 
C4 determine the potential V(a), we can use them to classify all possible cosmic types. 
All types of the universe in the theory of massive gravity are: 

(1) [Bounce]: If V(a) < for a G [ax, 00) and the equality holds at a = ax, the 
spacetime initially contracts from an infinite scale, and it eventually turns around at a 
finite scale ax, and then expands forever. 

(2) [Oscillation]: V(a) < for a G [a min , a max ] and the equality occurs at a = a min and 
a = CLmax- Thus, the spacetime oscillates between two finite scales. 

(3) [BB =3- BC\. V(a) < for a G (0, ax] and the equality holds at a = ax- The 
universe starts from a big bang (BB) and expands. Eventually it turns around at a = ax 
and contracts to a big crunch (BC). a T is the scale factor where the universe turns around 
from expansion to contraction. 

(4) [BB =^ 00 or 00 =>- BC]: V(a) < for any positive values of a. The spacetime 
starts from a big bang and expands forever, or the spacetime always contracts to a big 
crunch. 

III. EINSTEIN STATIC STATE SOLUTION 

Since we assume that the universe is dominated only by vacuum energy, p must be a 
positive constant, and thus A = > 0. The potential (Eq. (9)) can be re-expressed as 



V(a) 



2 

m 



(4c 3 + c 4 - 6 + 3A) a 3 + 



3(3 — 3c3 — C4) 



3a 



(4c 3 + c 4 - 6 + 3A) 




(10) 



5 



which means that V(a) = yields a cubic equation of a. An ES universe appears if there 
is a solution a = a s > which satisfies V(a s ) = and V'(a s ) = 0. At a s , both the speed 
of the cosmic expansion and acceleration are equal to zero and thus the universe can stay 
at this point in a long time if it is stable. Differentiating V(a) with respect to a, we have 

mP ( \ 

V\a) = -^2 (2(4c 3 + c 4 - 6 + 3A)a 3 + 3(3 - 3c 3 - c 4 )a 2 + (c 3 + c 4 ) J . (11) 

Combining V(a) = and V'(a) = 0, we find that an ES solution requires a relation 

between A and other two model parameters C3, c 4 to hold 

, 3 



A = A ± = ±2[1 + (C3 - l)c 3 + c 4 ]i - (c 3 - 1)[2 + c 3 (2c 3 - 1) + 3c 4 ] _ 

3(c 3 + c 4 ) 2 



They give two boundaries for obtaining an oscillating universe. Substituting Eq. (12) into 
the equation V(a) = or V'(a) = 0, one can find that the following static state solutions 

± 2c 3 + c 4 - 1 ± y/1 + ("I + c 3 )c 3 + c 4 

a s = a s = — — , 13 

3c 3 + c 4 - 3 

which is a double root of the equation V(a) = under the condition V'(a) = 0, and the 
third root is 

± c 3 + c 4 

Ot = a r = / = • (14) 

-1 + 2c 3 + c 4 ± 20 + (-1 + c 3 )c 3 + c 4 

If is positive, it corresponds to the radius where the universe turns around or bounces. 

Since the sign of 4c3 + c 4 — 6 + 3A, which is the coefficient of the a 3 term in the potential, 
plays a crucial role in determining the shape of the potential V(a), we will divide our 
following discussions into three different cases: 4c3 + c 4 — 6 + 3A > 0, 4c3 + c 4 — 6 + 3A < 
and 4c 3 + c 4 - 6 + 3A = 0. 



IV. THE CASE OF 4c 3 + c 4 - 6 + 3A > 

The condition 4c3 + c 4 — 6 + 3A > gives rise to a constraint on A, i.e., A > — |(4c3 + 
c 4 — 6). We find that, when A takes different values, the number of real roots for V(a) = 
is different. For example, A = A + allows the existence of three roots, where A + is defined 
in Eq. (12), but two of them are double, which corresponds to an unstable ES solution. 
In order to illustrate our results more clearly, we further divide our discussion into three 
subcases, i.e., A~ < A < A + , A < A~ and A > A + , respectively. 
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A. A" < A < A+ 



In this case, V(a) = yields a cubic equation of a, which has three real roots. Assuming 
a±, ci2 and a 3 are three real roots of this cubic equation, respectively, Eq. (10) can be 
expressed as 

2m 2 

Via) = (4c 3 + C4 — 6 + 3A)(a — aAia — 02) (a — a 3 ). (15) 

3a 

Since A > and (4c 3 + C4 — 6 + 3 A) > 0, the existence of three real roots requires that A 
must satisfy the condition 

Max jo, — i(4c 3 + c 4 — 6) j < A < A + , (16) 

if A~ < Max {0, -|(4c 3 + c 4 - 6)}, or 

A - < A < A+ , (17) 

if A~ > Max {0, — |(4c 3 + c 4 — 6)}. Since the cosmic scale factor must be larger than 
zero (a > 0), next, we will only consider the cases of positive roots. 



1. three positive roots 

Using a min , a max and ax to represent three positive roots of V(a) = and assuming 
< a min < a max < a T , we have 

Ira 2 

V{a) = — ^— (4c 3 + c 4 - 6 + 3 A) (a - a min )(a - a max )(a - a T ) 
3a 

2m 2 / 
= — 7T- (4c 3 + c 4 - 6 + 3A) a 3 - ia min + a max + a T )a 2 + 
3a \ 

Three positive roots mean that a min + a max + a T > 0, a min a max + a T a min + a max a T > 0, 
and a min a max a T > 0. Using these conditions and 4c 3 + c 4 — 6 + 3A > 0, and comparing 
Eq. (10) and Eq. (18), one can obtain three inequalities: 

3 - 3c 3 - c 4 < , c 4 + 2c 3 - 1 > , c 3 + c 4 > . (19) 



When the above inequalities are satisfied, A is negative or imaginary. Thus, A is required 
to only satisfy 



Max 



|-^(4c 3 + C4-6),o| <A< A + . 



(20) 



First, we consider the case of three different positive roots, which requires 
Max { — |(4c 3 + C4 — 6),0} < A < A + , and the A = A + case will be dealt with sep- 
arately. From this condition and Eq. (19), we find that three different positive roots 
demand 

3 



c 3 < - , c 3 + c 4 > 3 - 2c 3 
3 

or c 3 > - , c 3 + c 4 > . 




(21) 



FIG. 1: Potential V{a) under the conditions given in Eq. (21) and Max {-±(4c 3 + c 4 - 6), 0} < 
A < A + . An oscillating universe or a bouncing one is obtained. The constants are chosen as 
m = 1, c 3 = —5, c 4 = 20, and A = 2.1. The radii are a m j n = 0.64503, a max = 5.66024 and 
aT = 13.69470. The period of an oscillation is T = 5.52532. 

As shown in Fig. (1), V(a) < in a G [a min , a max ] with the equality holding at a = a min 
and a = a max , and a G [a^, 00) with V{aT) = 0. This means that the universe oscillates 
between a m i n and a max or bounces at ar- Thus, if the universe is in the region [a m j n , a max ] 
initially, it may undergo oscillation. After a number of oscillations, it may evolve to the 
bounce point ax through quantum tunneling but tunneling to the big bang singularity is 
not allowed. While, if the universe contracts initially from an infinite scale, it will turn 
around at a^ and then expand forever. So, the big bang singularity can be avoided in 
this case and a successful emergent universe is achieved. 
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FIG. 2: Potential V(a) with model parameters satisfying Eq. (21) and A = A + . There is an 
unstable static state universe and a bouncing one. The constants are chosen as m = 1, C3 = 1, 
c 4 = 1, and A = 0.47141. We obtain a min = 0.41421 and af = 3.41421. 

When A = A + in Eq. (20), a>T and a max coincide with each other, which corresponds 
to a double solution. We represent this double solution by af which is given in Eq. (13), 
and find that C3 and c 4 are required to satisfy Eq. (21), too. In Fig. (2), we plot the 
evolutionary curve of V(a) with A = A + . It is easy to see that af is an unstable ES 
solution. Therefore, the universe can oscillate between a m i n and af, and it can also 
evolve directly from a m j„ to 00 or evolve to 00 after some oscillations without the help of 
quantum tunneling. If the universe contracts initially from an infinite scale, it can turn 
around at af, or pass af and bounce at a m i n , then directly expand forever or do so after 
a number of oscillations between a m %n and af. 

In Fig. (3), we plot the phase diagram of spacetimes in (c 3 ,c 4 ) plane. Three positive 
roots restrict c 3 and c 4 to Region 1. 



2. two positive roots 

In this case, two of three real roots are positive and one of them is negative. We assume 
ai < < a2 < a%, and set ci2 = an and a% = 0^2 since a<i and 03 are two bouncing points. 
Then, we have 



aiaTiaT2 = (4c 3 + c 4 - 6 + 3A) < ° • (22) 
Of course, this condition corresponds to two different cases: three negative roots or two 
positive roots and one negative one. In order to distinguish these two cases, we further 
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FIG. 3: Phase diagram of spacetimes in (03,04) plane when A < A < A + . An oscillating 
universe is found in Region 1 and a bouncing one is obtained in Regions 1, 2, and 3. 



consider the signs of a\ + an + «T2 and a\an + an&T2 + a\ax2- When a\ + an + &T2 > 0, 
there is at least one positive root, which must correspond to the case of two positive 
roots and one negative one. When a\ + an + 0^2 < 0, three negative roots lead to 
a i a ri+^Ti a T2+ a i a T2 > 0, while, a\ < < an < Q-T2 implies that a 1 an+anaT2+ a iaT2 = 
«i(oti + 0LT2) + a T iO-T2 < —(an + «T2) 2 + a T idT2 = 
the conditions for two positive roots and a negative one are: 



(an + \(iT2) 2 - \a 2 T2 < 0. Thus, 



C3 + c 4 < , -(3 - 3c 3 - c 4 ) > , 
or c 3 + c 4 < , —(3 — 3c 3 — c 4 ) < , c 4 + 2c 3 — 1 < 



(23) 



Same as in the previous subsection, A is negative or imaginary when the above inequal- 
ities are satisfied, so, A is required to only obey 



Max I - ^(4c 3 + c 4 - 6), J < A < A + . 



(24) 



We first consider Max { — |(4c3 + C4 — 6), 0} < A < A + , which corresponds to the case 
of two different positive roots. Combing this condition with Eq. (23), we find that there 
exist two different positive roots and one negative one when 



c 3 < 1 , 3 - c 3 - 2^/3 - 3c 3 + c\ < c 3 + c 4 < , 
^ < c 3 < 2 , 3 - 2c 3 < c 3 + c 4 < , 



or c 3 > 2 , 3 — c 3 — 2 \l 3 — 3c 3 + cj < c 3 + c 4 < . 



(25) 
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FIG. 4: The potential V{a) for a 55 =>• BC universe or a bouncing one with Max {— g(4c 3 + 
C4 — 6),0} < A < A + and model parameters in Region 2 of Fig. (3). The constants are chosen 
as m = 1, c 3 = 3, C4 = —5, and A = 0.1. The radii are clti = 1.15556 and ax2 = 1-86572. 

Above conditions correspond to Region 2 in Fig. (3) . 

Fig. (4) shows the evolution of the potential V(a) with the model parameters in Region 
2 of Fig. (3). From this figure, we find that a BB =>- BC universe or a bouncing one 
is obtained since V < in a G (0, an] and [ax2,oo) with V = occurring at a = axi 
and a = 0^2 • Thus, if the universe initiates from a big bang, it can expand to an- It 
then turns over at an and ends with a big crunch. In addition, it is also possible that 
the universe quantum tunnels to a^2 directly from a^i and then expands forever. If the 
universe contracts initially from infinity, a bounce will occur at ot2- 

When A = A + , a^i and a<r 2 coincide with each other and a double root a+ given in 
Eq. (13) is obtained. c 3 and c 4 are required to satisfy Eq. (25), too. As shown in Fig (5), 
the ES solution a+ is unstable. So, if the universe initiates from big bang, it will expand 
to an unstable ES universe and then turn over or expand forever. If the universe contracts 
from infinity initially, it can bounce at a+ or end with a big crunch. 

3. one positive root 

For this case, only one of three real roots is positive and two of them are negative. 
Assuming that a 1; a 2 < and a 3 = a T > 0, one has 




> 0. 



(26) 



(4c 3 + c 4 - 6 + 3A) 
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FIG. 5: Potential V{a) with c 3 and c 4 in Region 2 of Fig. (3), and A = A + . A big bang universe 
and an unstable ES solution are obtained. The constants are chosen as m = 1, c 3 = 3, c 4 = —6.3 
and A = 0.15217. We find a+ = 1.54447. 

Since this condition can correspond to either three positive roots or one positive root and 
two negative ones, we have to consider other conditions coming from a\ + a 2 + «t and 
aia 2 + a 2 aT + aidr- If a\ + + ar < 0, there is at least one negative root regardless of 
the value of a 4 a 2 + a 2 a^ + cl^t- Thus, 



are sufficient for having one positive root and two negative ones. If a\ + a 2 + clt > 0, the 
condition from ai<2 2 + a2«T + clicit should be added. When a±, a 2 , ax > (three positive 
roots), aia 2 + a 2 ay + a\ar > 0. However, if a\ < a 2 < < a^, then aia 2 + a 2 a^ + cliclt = 
a T (ai + a 2 ) + a x a 2 < -(ai + a 2 ) 2 + a x a 2 = -(a x - \a 2 ) 2 - |a 2 < 0. So, 



give the conditions in obtaining one positive root and two negative ones. Since A < if 
Eq. (27) or Eq. (28) is satisfied, we only consider 



Here, A = A + is discarded because in this case there is no double solution. That is, there 
is no unstable ES solution. Combining Eqs. (27) and (29), or Eqs. (28) and (29), one can 
find the conditions for the existence of one positive root and two negative ones 



c 3 + c 4 > , -(3 - 3c 3 — C4) < , 



(27) 



c 3 + c 4 > , -(3 - 3c 3 - c 4 ) > , c 4 + 2c 3 - 1 < , 



(28) 




(29) 



c 3 < - , < c 3 + c 4 < 3 - 2c 3 , 



(30) 
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which is shown as Region 3 in Fig. (3). Fig. (6) shows the evolutionary curve of the 
potential V(a) with model parameters in Region 3 in Fig. (3). Apparently, a bouncing 
universe is obtained. 
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FIG. 6: Potential V(a) with Max {— g(4c3 + C4 — 6),0} < A < A + and model parameters in 
Region 3 of Fig. (3). A bouncing universe is obtained. The constants are chosen as m = 1, 
C3 = —3, C4 = 6, and A = 10. The turning radius at a bounce is clt = 0.408248. 



4- no positive root 

Three real roots are all negative. Fig. (7) shows the evolutionary curve of V(a) in this 
case. One can see that the potential is always negative and the cosmic evolution type is 
BB =>- 00 or 00 =>- BC. 




FIG. 7: Potential V(a) for the case of no positive root. A BB => 00 or 00 => BC cosmic type 
is obtained. The constants are chosen as m = 1, C3 = 1.5, C4 = —1.7, and A = 0.58. 
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B. A > A+ 



In this case, A must satisfy 



A > Max < A + , --(4c 3 + c 4 - 6), 



(31) 



Under this condition, there is only one real root, which can be positive or negative, and 
other two roots are a conjugate imaginary pair. In Fig. (8), we show all possible cosmic 
types in (c 3 , c 4 ) plane. 



FIG. 8: Phase diagram of spacetimes in (03,04) plane with A > A + . A bouncing universe is 
found in Region 1, while, a BB BC one is obtained in Region 2. 

1. one positive root 

Assuming that a\ and 02 are two conjugate imaginary roots and 03 is the only positive 
one, and setting a 3 = a^, one has 



Combining this inequality with Eq. (31), we get the condition for one positive root and 
two conjugate imaginary ones 





> 0. 



(32) 



(4c 3 + c 4 - 6 + 3A) 



c 3 + c 4 > 0, 



(33) 



14 



which corresponds to Region 1 of Fig. (8). Fig. (9) shows the evolution of the potential 
V(a) with model parameters in Region 1 of Fig. (8). We find that a bouncing universe is 
obtained. 




FIG. 9: Potential V(a) for a bouncing universe with model parameters in Region 1 of Fig. (8) 
and A > A + . The constants are chosen as m = 1, c 3 = 1.5, c 4 = 2, and A = 2. The turning 
radius at a bounce is ax = 0.3823023. 



2. no positive root 

In this case, the only real root is negative, which gives 

(c 3 + C 4 ) 

aia2fl3 = (4c 3 + c 4 - 6 + 3A) < ° ' (34) 
From the above inequality and the condition given in Eq. (31), one has 

c 3 ^l, -(C 3 -1) 2 <C3 + C4<0, (35) 

which correspond to Region 2 in Fig. (8). Fig. (10) shows that the potential is always 
negative and the cosmic type is BB 00 or 00 =>■ BC. 

C. A < A- 

It is easy to see that here A must satisfy 

Max jo, — -^(4c 3 + c 4 — 6) j < A < A". (36) 
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FIG. 10: Potential V(a) for a BB =;> oo or oo => BC universe with A > A + and model 
parameters in Region 2 of Fig. (8). The constants are chosen as m = 1, c 3 = 1.5, C4 = —1.6, 
and A = 2. 

One can show that there is only one real root and it is negative. Other two roots are a 
conjugate imaginary pair. Model parameters c 3 and C4 are restricted in the region: 



1 < c 3 < 2 , -(c 3 - I) 2 < c 3 + c 4 < 3 - c 3 - 2a 3 - 3c 3 + c 



(37) 



The evolution of the potential is shown in Fig. (11). One can see that the cosmic type is 
BB =>• 00 or 00 =3- BC. 

In Tab. (I), we sum up the results obtained in this section. 



0.5 1.0 1.5 2.0 2.5 a 




FIG. 11: Potential V(a) under the condition {0, — |(4c 3 + C4 — 6)} < A < A~ and model 
parameters satisfying Eq. (37). A BB =4> 00 or 00 => BC universe is obtained. The constants 
are chosen as m = 1, c 3 = 1.5, C4 = —1.74, and A = 0.6. 
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TABLE I: Summary of the cosmic type in the 4c 3 + C4 — 6 + 3A > case 



A 


C3, c 4 


Cosmic Type 


Max {-|(4c 3 + c 4 + 6), 0} 
< A< A+ 


C3 < |, c 3 + c 4 > 3 - 2c 3 
or c 3 > | , c 3 + c 4 > 


Oscillation 
or Bounce 


c 3 < 1, 3 - c 3 - 2^3 - 3c 3 + c 2 , < c 3 + c 4 < 
| < c 3 < 2, 3 - 2c 3 < c 3 + c 4 < 
or c 3 > 2, 3 — c 3 — 2^3 — 3c 3 + c\ < c 3 + c 4 < 


BB BC 
or Bounce 


c 3 < §, < c 3 + c 4 < 3 - 2c 3 


Bounce 


A = A+ 


C3 < |, c 3 + c 4 > 3 - 2c 3 
or c 3 > | , c 3 + c 4 > 


Bounce 
Unstable ES 


c 3 < 1, 3 — C3 — 2^ 3 - 3c 3 + c| < c 3 + c 4 < 
| < c 3 < 2, 3 - 2c 3 < c 3 + c 4 < 
or c 3 > 2, 3 - c 3 - 2^3 - 3c 3 + c| < c 3 + c 4 < 


BB => 00 
=> BC 
00 => BC 
or Bounce 
Unstable ES 


A > 

Max{A + ,-i(4c 3 + c 4 -6),0} 


c 3 + c 4 > 


Bounce 


c 3 / 1, -(C3 - l) 2 < c 3 + c 4 < 


BB => 00 
or 00 =>■ BC 


Max {0,-±(4c 3 + c 4 - 6)} 
< A < A" 


K c 3 < 2, 


BB => 00 
or 00 =>• BC 


-(c 3 - l) 2 < c 3 + c 4 < 3 - c 3 - 2^3 - 3c 3 + c{ 



V. THE CASE OF 3A + 4c 3 + c 4 - 6 < 

As in the preceding section, we divide our discussion into three different subcases: 
A~ < A < A+, A < A" and A > A + , respectively. 
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A. A" < A < A+ 



In this case, the equation V(a) = has three real roots and we assume them to be a\, 
a 2 and a 3 , respectively. Since a positive A is considered, A must satisfy 

< A < --(4c 3 + c 4 - 6) (38) 

and 

A~ < A < A + . (39) 

1. three positive roots 

Letting a\ = a T , a 2 = a min , and a 3 = a max , and assuming < a T < a min < a max , we 
can re-express Eq. (10) as 

2m 2 

v ( a ) = — ^— (4c3 + c 4 - 6 + 3A)(a - a T )(a - a min )(a - a max ) 

2m 2 ( 
= — ^~( 4c3 + c 4 - 6 + 3A) ( a 3 - (a T + a min + a max )a 2 + 

u 'iriax I • 

(40) 



Three positive roots imply that a T + a min + a max > 0, a T a min + a min a max + a T a max > 
and aTdminO'max > 0. Comparing Eq. (10) and Eq. (40), one has 

3 - 3c 3 - c 4 > , c 4 + 2c 3 - 1 < , c 3 + c 4 < 0. (41) 

We first study the case of A ^ A + and A ^ A - , which corresponds to three different 
positive roots. Combining Eqs. (38, 39) and Eq. (41), we obtain that three different 
positive roots require 

3 

c 3 < 1 , --(c 3 - l) 2 < c 3 + c 4 < , 
2 < c 3 < 3 , -(c 3 - l) 2 < c 3 + c 4 < 3 - 2c 3 , 
or c 3 > 3 , -(c 3 - l) 2 < c 3 + c 4 < 6 - 3c 3 , (42) 

which give Region 1 in Fig. (12). The cosmic type is BB =>- BC or oscillation, as can 
be seen from Fig. (13). Thus, if the universe starts from a big bang, it can expand to 
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a^, then turn over at a? and end with a big crunch. If the universe is in the region 
a G [dmini o-max ] initially, it may undergo an oscillation. After some oscillations, it may 
quantum mechanically tunnel to a<r and end with a big crunch singularity. Therefore, the 
classical singularity still exists. 
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FIG. 12: Phase diagram of spacetimes in (03,04) plane under the condition of A~ < A < A + . 
An oscillating universe is found in Regions 1 and 2, and a BB =>- BC universe is obtained in 
Regions 1 and 3. 
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FIG. 13: Potential V(a) for a BB =>• BC universe or an oscillating one with model parameters 
in Region 1 of Fig. (12). The constants are chosen as m = 1, C3 = 2.5, C4 = —4.6, and A = 0.196. 
The radii are aj- = 1.53548, o m m = 

6.86654 and 

a max — 16.598. The period of an oscillation is 

T = 107.478. 

When A = A + and A satisfies < A < — |(4c3 + C4 — 6), a? and a m i n coincide with 
each other and forms a double positive root a+ defined in Eq. (13). Fig. (14) shows that 
a+ is an unstable ES solution. The cosmic type is BB =>- BC. But, the universe can turn 
over at af, or a max . From Eq. (40) and < A = A + < — |(4c3 + C4 — 6), we obtain the 
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conditions on c 3 and C4 for an unstable Einstein static state solution: 

c 3 < 1 , --(c 3 - l) 2 < c 3 + c 4 < 3 - c 3 - 2^3-3c 3 + c§ , 
or c 3 >2, -(c 3 - l) 2 < c 3 + c 4 < 3 - c 3 - 2v/3 - 3c 3 + c\ . (43) 
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FIG. 14: The potential V(a) under the condition A = A + . An unstable ES solution is obtained. 
The constants are chosen as m = 1, c 3 = —1, C4 = —1, and A = 0.804738. The radii are 
a+ = 0.369398 and a max = 1.70711. 

If < A = A~ < — |(4c 3 + c 4 — 6), there is also a double solution a~, which is the 
coincidence of a m i n and a max , and, as shown in Fig. (15), it is a stable ES solution. This 
stable ES solution requires that c 3 and C4 satisfy 

2<c 3 <3, -(c 3 -l) 2 <c 3 + c 4 <3-2c 3 , 
or c 3 >3, -(c 3 -l) 2 <c 3 + c 4 <-Jj(c 3 -l) 2 . (44) 

Since V(a) < in a G (0, a^] and a = a~ , the cosmic type is BB BC if the scale 
factor is less than ar initially, or the universe stays at a~ . However, the universe can not 
stay at this stable ES past-eternally since quantum tunneling may drive it into the region 
a G (0, a?}. Thus, the big bang singularity can not be avoided although there is a stable 
ES solution. As a result, in massive gravity the existence of a stable ES solution can not 
successfully resolve the big bang singularity problem. 

In addition, there is also a possibility such that A = A + = A - , which means that ax, 
a m i n and a max merges to form a triple root a^. Using 

0< A = A + = A~ < -i(4c 3 + c 4 -6) , (45) 

o 
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and Eq. (41), one can find that the conditions for a triple root are C3 > 2, C4 = — 1 + C3 — C3 
and A = 3( - e3 1 _ 1 ) ■ From Fig. (16), we can see that the potential V(a) < when a £ (0, ar], 
which means that the cosmic type is BB =>- BC . 




FIG. 15: Potential V(a) under the condition A = A .A BB =>■ BC universe or a stable ES one 
is obtained. The constants are chosen as m = 1, C3 = 2.5, C4 = —4.6, and A = 0.195299. The 
radii are a T = 1.52772 and a~ = 9.87298. 




FIG. 16: Potential V(a) under the condition A = A + = A . A BB =>• BC universe is obtained. 
The constants are chosen as m = 1, C3 = 3, C4 = —7, and A = 0.166667. The radius is ay = 2. 



2. two positive roots 



We assume a\ < 0, and 02, 03 > 0, and set a-i = a m i n and 02 = a max , so that 



aiammflmox = Ta i r I QA\ < ' ( 46 ) 

(4c 3 + c 4 — 6 + 3A) 



21 



An analysis similar to that in the previous section leads to the conditions for two positive 
roots and one negative root as follows: 

c 3 < 2 , < c 3 + c 4 < 6 - 3c 3 , (47) 

which determine Region 2 of Fig. (12). We find that A = A^ 1 is forbidden since A + > 
— |(4c 3 + C4 — 6) and A~ < when c 3 and C4 satisfy Eq. (47). This implies that A is only 
required to satisfy Eq. (38). Fig. (17) shows the evolution of the potential V(a) with the 
model parameters in Region 2 of Fig. (12). From Fig. (17), we find that an oscillating 
universe is achieved since V < in a e [a min , a max ] with V = occurring at a = a min and 

&max • 
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Potential V{a) for an oscillating universe with model parameters in Region 2 of 
The constants are chosen as m = 1, c 3 = 1, C4 = —0.5, and A = 0.5. The radii 
0.27255 and a max = 2.10074. 



FIG. 17: 
Fig. (12). 

are drain — 



3. one positive root 

Assuming that 01,02 < and a 3 = a<r > and using the analysis similar to that in 
the previous section, we get that the conditions for one positive root 

l<c 3 <2, 3 - c 3 - 2^3 - 3c 3 + c§ < c 3 + c 4 < , 
or 2 < c 3 < 3 , 3 - 2c 3 < c 3 + c 4 < 6 - 3c 3 . (48) 

Region 3 of Fig. (12) shows the allowed values of c 3 and c 4 for only one positive root. 
Fig. (18) shows the evolution of the potential V(a). It is easy to see that a BB — > BC 
universe is realized. 
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FIG. 18: Potential V(a) for a BB =>- BC universe with model parameter in Region 3 of Fig. (12). 
The constants are chosen as m = 1, c 3 = 2.5, C4 = —4.25, and A = 0.06. and the turning radius 
is a T = 1.07118. 

4- no positive root 

Since a±, a 2 , a 3 < 0, and 4c 3 + c 4 — 6 + 3 A < 0, one has three inequalities: 

3 - 3c 3 - c 4 < 0, c 4 + 2c 3 - 1 < 0, c 3 + c 4 > 0. (49) 

When c 3 and c 4 satisfy the above inequalities, A~ < and A + < 0. Therefore, there is no 
allowed positive value for A. This means that this is not a physically meaningful case. 



B. A > A+ 

This case corresponds to only one real root a 3 . Other two roots a 4 and a 2 are a 
conjugate imaginary pair. Now A must satisfy 

Max {0,A + } < A< -i(4c 3 + c 4 -6) (50) 

If this real root a 3 is negative, then 

010203 = Ta , a i QA'i < °- 51 

(4c 3 + c 4 - 6 + 3A) 

Considering the condition of A (Eq. (50)), we find that there is no solution for Eq. (51). 
Thus, a 3 must be a positive one and we set a 3 = a^, which means 

a ' a ^= (4c3^-6 ) + 3A) > °- < 52 > 
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FIG. 19: Phase diagram of spacetimes in (c 3 , C4) plane when A > A + . A BB =>• BC universe is 
found in the gray region. 

Combining Eq. (50) and Eq. (52), we obtain that c 3 and c 4 must obey 

c 3 <l, -(c 3 - l) 2 < c 3 + c 4 < 3 - c 3 - 2^3 - 3c 3 + cj , 
or c 3 >2, -(c 3 - l) 2 < c 3 + c 4 < 3 - c 3 - 2^3 - 3c 3 + c§ , (53) 

in order to get one positive root and two imaginary roots which are conjugate to each 
other. In Fig. (19), we show the allowed values of c 3 and C4 for a positive root. Fig. (20) 
shows the evolution of the potential V(a) with the model parameters in the gray regions 
of Fig. (19), and a BB BC universe is obtained. 




FIG. 20: Potential V{a) with model parameters in the gray regions of Fig. (19) and A > A + . 
We obtain a BB =4> BC universe. The constants are chosen as m = 1, c 3 = 3, C4 = —6.75, and 
A = 0.2. The radius is a T = 12.3248. 
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C. A < A~ 



This case requires 

< A < Mini A - , — ^(4c 3 + c 4 — 6)1 . (54) 

Under this condition Eq. (51) can't be satisfied, but Eq. (52) can. Thus, there is also 
one positive root a 3 = a? and two conjugate imaginary ones a\ and a 2 . Combining the 
condition on A given in Eq. (54) and Eq. (52), we get 

c 3 > 1 , -(ca - l) 2 < c 3 + c 4 < --(ca - l) 2 . (55) 

The allowed region of c 3 and c 4 is shown in Fig. (21) and the evolutionary curve of the 
potential is given in Fig. (22). We find that the potential V(a) < when < a < ax and 
the cosmic type is BB =>• BC. 

Tab. (II) sums up the results of this section. 

' ' ^ ' ' ' ' ^ ' ' ' ' 4 ' ' ' ' 5 c 3 




FIG. 21: Phase diagram of spacetimes in (03,04) plane under the condition A < A . A BB 
BC universe is found in the gray region. 



VI. THE CASE OF A = -i(4c 3 + c 4 -6) 



In this case Eq. (9) becomes 



V(a) = ™ 
a 



(3 - 3c 3 - c 4 )a 2 - (c 4 + 2c 3 - l)a + -(c 3 + c 4 ) 



(56) 
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FIG. 22: Potential V(a) with model parameters in the gray region of Fig. (21) and A < A . 
We obtain a BB =4> BC universe. The constants are chosen as m = 1, c 3 = 2, c 4 = —2.85, and 
A = 0.15. The radius is a T = 1.22118 



Apparently, V(a) = has two roots: 



3 — 6c 3 — 3c 4 — a/3 (3 + 6c 4 — 4c 3 c 4 — cf) 
1 —6(3 — 3c 3 — c 4 ) 

_ 3 - 6c 3 - 3c 4 + a/3 (3 + 6c 4 - 4c 3 c 4 - c 4 ) 
a2 ~ -6(3 - 3c 3 - c 4 ) (57) 

Now we divide our discussion into two cases, i.e., — (3— 3c 3 — c 4 ) > and — (3— 3c 3 — c 4 ) < 0. 
A. -(3-3c 3 -c 4 ) >0 




C3 



FIG. 23: Phase diagram of spacetimes in (c 3 , c 4 ) plane under the condition A = — 3 (4c 3 + c 4 — 6) 
and —(3 — 3c 3 — c 4 ) > 0. An oscillating universe is found in Region 1, and a BB BC one is 
obtained in Region 2. 
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TABLE II: Summary of the cosmic type in the 4c 3 + C4 — 6 + 3A < case 



A 


c 3 , c 4 


Cosmic Type 


Max {A",0} < A < 
Mm{-±(4c 3 + c 4 + 6),A+} 


c 3 < 1, -|(c 3 - l) 2 < c 3 + c 4 < 
2 < c 3 < 3, -(c 3 - l) 2 < c 3 + c 4 < 3 - 2c 3 
or c 3 > 3, -(c 3 - l) 2 < c 3 + c 4 < 6 - 3c 3 


BB => BC 
or Oscillation 


c 3 < 2, < c 3 + C4 < 6 — 3c 3 


Oscillation 


1 < c 3 < 2, 3 - c 3 - 2^/3 - 3c 3 + c\ < c 3 + c 4 < 
or 2 < c 3 < 3, 3 — 2c 3 < c 3 + C4 < 6 — 3c 3 


BB BC 


< A = A+ < 
-s(4c3 + c 4 + 6) 


c 3 < 1, -| (c 3 - l) 2 < c 3 + c 4 

< 3 - c 3 - 2^/3 - 3c 3 + 

or c 3 > 2, -(c 3 - l) 2 < c 3 + c 4 

< 3 - c 3 - 2^3 - 3c 3 + c 2 3 


BB => BC 
Unstable ES 


< A = A" < 
-|(4c 3 + c 4 - 6) 


2 < c 3 < 3, -(c 3 - l) 2 < c 3 + c 4 < 3 - 2c 3 
or c 3 > 3, -(c 3 - l) 2 < c 3 + c 4 < -|(c3 - I) 2 


BB => BC 
Stable ES 


A " A ~ " A+ " 3( C 3-1) 


c 3 > 2, c 3 + c 4 = -(c 3 - l) 2 


BB BC 


Max {0, A+} < A < 
-±(4c 3 + c 4 -6) 


c 3 < 1, -(c 3 - l) 2 < c 3 + c 4 

< 3- 03-2^3-303 + ^ 
or c 3 > 2, -(c 3 - l) 2 < c 3 + c 4 

< 3 - c 3 - 2^3 - 3c 3 + c 2 3 


BB => BC 


< A < 
Min{A~, -|(4c 3 + c 4 - 6)} 


c 3 > 1, -(c 3 - l) 2 < c 3 + c 4 < -|(c 3 - l) 2 


BB => BC 



If two roots are all positive, we set 02 = a m i n and a 3 = a max , and find that c 3 and 
C4 must satisfy A = (3 + c 4 (6 — 4c 3 ) — c 2 ) > 0, c 3 + c 4 > 0, —(3 — 3c 3 — c 4 ) > 0, 
A = -|(4c 3 + c 4 - 6) > 0, and -(c 4 + 2c 3 - 1) < 0. These lead to 

3 

c 3 < - , 3 - 2c 3 < c 3 + c 4 < 6 - 3c 3 , 

3 , v 

or - < c 3 < 2 , < c 3 + c 4 < 6 - 3c 3 , (58) 
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which is represented as Region 1 in Fig. (23) where the phase diagram of spacetimes 
in (03, C4) plane is shown. Fig. (24) displays the evolution of the potential with model 
parameters in Region 1 of Fig. (23). We can see that V(a) < in a e [a m j n ,a max ] with 
the equality holding at a = a m i n and a = a max . Thus, an oscillating universe is obtained. 

If a\ = a 2 , then C4 = 3 — 2c 3 ± 2 a/3 + C3 (03 — 3) , which is outside Region 1 of Fig. (23). 
Thus, a stable static Einstein universe can't be obtained in this case. 




2.0 2.5 a 



FIG. 24: Potential V(a) with model parameters in Region 1 of Fig. (23). An oscillating universe 
is obtained. The constants are chosen as m = 1, C3 = 1, C4 = 1, and A = |. The radii are 
o-min = 0.42265 and a max = 1.57735. 



Now we consider the case of a\ < and a<i = clt > 0. We find that there is a BB 
cosmic evolution type as shown in Fig. (25), and c 3 and c 4 satisfy 



BC 



- < c 3 < 2 
2 



3 - 2c 3 < c 3 + c 4 < 



or 2 < C3 < 3 , 3 — 2C3 < C3 + C4 < 6 — 3c 3 



(59) 



which correspond to Region 2 in Fig. (23). 

If both a\ and a 2 are negative, then V(a) > in a G (0, 00). So this is not a case of 
physical significance. 
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FIG. 25: Potential V(a) with model parameters in Region 2 of Fig. (23). A BB => BC universe 
is obtained. The constants are chosen as m = 1, C3 = 2.5, C4 = —4.25, and A = 0.083333. and 
the turning radius is ax = 1.10728. 
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FIG. 26: Phase diagram of spacetimes in (03, C4) plane under the conditions A = — |(4c3 + C4 — 6) 
and —(3 — 3c3 — C4) < 0. A bouncing universe is found in Regions 1 and 2, and a BB =>■ BC 
universe is obtained in Region 1. 

B. -(3 - 3c 3 - c 4 ) < 

Since a\ > and 02 > require A = |(3 + c 4 (6 — 4c 3 ) — cf) > 0, c 3 + c 4 < 0, and 
— (04 + 2C3 — 1) > 0, the conditions for two positive roots are 

c 3 < 1 , 3 - c 3 - 2^3 - 3c 3 + c§ < c 3 + c 4 < , 
2<c 3 <3, 3-C3- 2^3-303 + 01 < c 3 + c 4 < 3 - 2c 3 , 
or c 3 > 3 , 3 — c 3 — 2yJ 3 — 3c 3 + c| < c 3 + c 4 < 6 — 3c 3 . (60) 

In Region 1 of Fig. (26), we show the allowed values C3 and C4 for two positive roots. 
Setting ai = a T2 and a 2 = ari since a x > a 2 , we show the evolution of the potential V(a) 
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in Fig. (27) with the model parameters in Region 1 of Fig. (26). From this figure, we 
see that a BB =>• BC universe or a bouncing one is obtained since V < in a G (0, an] 
and [0^2, oo) with V = occurring at a = an and a = aj- 2 . Another possibility is that 
the universe starts at the big bang singularity and, rather than bounce back, it quantum 
mechanically tunnels to an when it evolves to aT2 and then expands forever. 




FIG. 27: Potential V(a) with model parameters in Region 1 of Fig. (26). A BB => BC universe 
or a bouncing one is found. The constants are chosen as m = 1, c 3 = 2.5, c 4 = —4.55, and 
A = 0.183333. The radii are a T i = 1.42774 and a T 2 = 9.57226. 



If c 4 = 3 — 2c 3 ±2a/3 + c 3 (c 3 — 3), it is easy to see a x = a 2 , which means that V(a) = 
has a double solution a. 



-3 + 2c 3 + 2^3 + c 3 (c 3 - 3) 

O's = s • 

3(c 3 - 2) 

We find that a s is an unstable ES universe and it exists under the conditions 



(61) 



c 3 < 1 , c 3 + c 4 = 3 - c 3 - 2a/3 + c 3 (c 3 - 3) , 
or c 3 > 2 , c 3 + c 4 = 3 - c 3 - 2 a/3 + c 3 (c 3 - 3) . 



(62) 



Thus, as shown in Fig (28), if the universe initiates from big bang, it will expand to an 
unstable ES universe and then turn over or expand forever. 

Now we consider the case of a 2 < and a± = (It > 0. In this bouncing; universe 

is obtained, as shown in Fig. (29), if c 3 and c 4 satisfy 



c 3 < - , < c 3 + c 4 < 3 - 2c 3 , 



(63) 



which corresponds to Region 2 of Fig. (26). 
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FIG. 28: Potential V(a) with model parameters satisfying Eq. (63). An unstable Einstein static 
state universe is obtained. The constants are chosen as m = 1, C3 = 3, C4 = —4.64575, and 
A = 0.15470. a s = 1.57735. 



FIG. 29: Potential V{a) with model parameters in Region 1 of Fig. (26). A bouncing universe 
is found. The constants are chosen as m = 1, C3 = 1, C4 = —0.5, and A = 0.83333. The turning 
radius at a bounce is ax = 0.263763. 

If both a>i and a 2 are not positive, we find that the potential is always negative as 
shown in Fig. (30) and the cosmic type is BB =>- 00 or 00 =r- BC. 
The results of this section are summed up in Tab. (III). 

VII. CONCLUSIONS 

Massive gravity is a modification of general relativity. It has been spurring an increasing 
deal of interest recently, since it can explain the present accelerated cosmic expansion 
without the need of dark energy. In this paper, using a method in which the scale factor a 
changes as a particle in a "potential" , we analyze all possible cosmic evolutions in a ghost- 
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FIG. 30: Potential V{a) for the case of without positive root. A BB =4> oo or oo =>• BC one is 
obtained. The constants are chosen as m = 1, c 3 = 1.2, c 4 = —1.22, and A = 0.806667. 



TABLE III: Summary of the cosmic type in the 4c 3 + C4 — 6 + 3A = case 



c 3 , c 4 



Cosmic Type 



C3 < |, 3 - 2c 3 < c 3 + c 4 < 6 - 3c 3 
or § < c 3 < 2, < c 3 + c 4 < 6 - 3c 3 



Oscillation 



I < c 3 < 2, 3 - 2c 3 < c 3 + c 4 < 



or 2 < c 3 < 3, 3 — 2c 3 < c 3 + c 4 < 6 — 3c 3 



BB => BC 



c 3 < 1, 3 - c 3 - 2 1/3 - 3c 3 + c| < c 3 + c 4 < 



2 < c 3 < 3, 3 - c 3 - 2^3 - 3c 3 + cj < c 3 + c 4 < 3 - 2c 3 



SB => BC 
or Bounce 



or c 3 > 3, 3 — c 3 — 2 a/3 — 3c 3 + c| < c 3 + c 4 < 6 — 3c 3 



c 3 < 1, c 3 + c 4 = 3 - c 3 - 2i/3 + c 3 (c 3 - 3) 



or 



c 3 > 2, c 3 + c 4 = 3 - c 3 - 2^3 + c 3 (c 3 - 3) 



BB => BC 
BB^ 00 
00 =>■ BC 
or Bounce 
Unstable ES 



c 3 < |, 0<c 3 + c 4 <3-2c 3 , 



Bounce 



free massive gravity theory. A spatially flat universe is considered in our discussion and we 
assume that the vacuum energy is the only energy component. The results are summed 
up in Tabs. (I, II, III). We find that there may exist an oscillating universe between a m i n 
and a max or a bouncing one at if model parameters are in some specific regions. If 
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the cosmic scale factor is in the region [a min , a max ] initially, the universe may undergo an 
oscillation. After a number of oscillations, it may evolve to the bounce point through 
quantum tunneling. While, if the universe contracts initially from an infinite scale, it 
can turn around at ar and then expand forever. Thus, the big bang singularity problem 
can be resolved successfully. Remarkably, although we do have a stable ES solution in 
some circumstances, the universe can not stay at this stable state past-eternally since it 
is allowed to quantum mechanically tunnel to a big-bang-to-big-crunch cosmic evolution 
type and end with a big crunch. Thus, the existence of a stable ES universe can not 
successfully resolve the big bang singularity in the massive gravity. This feature is related 
to the behavior of V(a) — > — oo as a — > 0, which is a result of C3 + C4 < in the massive 
gravity we consider in the present paper, when there exists a stable ES universe with a 
finite a s . Let us note however that both in the Horava-Lifshitz gravity [18] and the DGP 
braneworld scenario [19], there exist stable ES universes where V(a) — > 00 as a — > 0, so 
a quantum tunneling to the big bang singularity is not allowed ( refer to Fig. (4) in [18] 
and Fig. (5) in [19]) and as a result the existence of a stable ES universe can avoid the 
big bang singularity in these theories. Therefore, whether the existence of a stable ES 
universe can resolve the big bang singularity or not is a peculiar feature of the theory of 
gravity itself and it is in fact determined by the behavior of the leading term in V(a) as 
a approaches zero, i.e., the big bang singularity. 
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